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Abstract 

Gamma-ray telescopes have reported some surprising observations of multi- 
TeV photons from distant active galactic nuclei (AGN), which show no sig- 
nificant attenuation due to pair production on either the extragalactic back- 
ground light (EBL), or the photons near the source. We suggest a new 
interpretation of these observations, which is consistent with both the EBL 
calculations and the AGN models. Cosmic rays with energies below 50 EeV, 
produced by AGN, can cross cosmological distances, interact with EBL rela- 
tively close to Earth, and generate the secondary photons observed by 7-ray 
telescopes. We calculate the spectrum of the secondary photons and find that 
it agrees with the 7-ray data. The delays in the proton arrival times can ex- 
plain the orphan flares, the lack of time correlations, and the mismatch of the 
variability time scales inferred from the multiwavelength observations. The 
7-ray data are consistent with the detection of the secondary photons, which 
has important ramifications for gamma-ray astronomy, cosmic ray physics, 
EBL, and the intergalactic magnetic fields (IGMF). 
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The recent observations of multi-TeV photons from distant blazars [H, 0, 0] 
are surprising because their spectra do not show the expected suppressions 
due to the interactions with extragalactic background light (EBL). The ob- 
servations have led to a debate as to whether one should reconsider the 
calculations of EBL 0, 0, @, 0]. However, even if EBL is thinner and softer, 
the observations of very high energy (VHE) photons are difficult to recon- 
cile with the expected absorption at or near the source. For example, the 
intrinsic photon densities around M87 and 3C279 can be inferred from ob- 
servations, and these densities appear inconsistent with the unattenuated 
emission of VHE photons 0, @| . Among other proposed explanations, it was 
suggested that photons may convert into some hypothetical axion-like par- 
ticles that convert back into photons in the galactic magnetic fields @,llq| 



or that Lorentz invariance violation may be the reason for the lack of ab- 



sorption The orphan flashes, the lack of timing correlations between 



VHE and optical observations 12] and the mismatch of time scales for X-ray 



and VHE variabilities [131 ] . as well as some unusual energy-dependent time 
delays [l4j, pose additional questions. 

These puzzling inconsistencies suggest an alternative interpretation of the 
observations of distantly] blazars. Active galactic nuclei (AGN) are believed to 
produce ultrahigh-energy cosmic rays. For energies below Greisen-Zatsepin- 
Kuzmin (GZK) cutoff (i.e., below 50 EeV), the cosmic rays can propagate 
cosmological distances without interacting. The deflection of their trajec- 
tories depends on the intergalactic magnetic fields (IGMF), which can be 
very small along the line of sight (only the upper limits can be derived from 
observations). If the IGMF small as predicted by the constrained sim- 

ulations [l7], the trajectories of high-energy protons are close to the line of 
sight until the protons enter our galaxy, where the magnetic fields are much 



strongeip. For cosmic rays with energies (4 — 8) x 10 eV, a correlation 



of the arrival directions with BL Lacertae objects has been reported [20 
For lower energies, the deflection of protons in the galactic magnetic fields 
makes the association with an extragalactic source impossible. However, 
the secondary photons produced by lower-energy cosmic rays point back to 



1 For nearby AGN, the primary photons are expected to dominate the signal, and the 
EBL absorption features should be seen Til, 1(| ■ We concentrate on the AGN that are far 
enough for the secondary photons to be important. 

2 These results depend on a number of assumptions, 



ultimately be the best probe of IGMF [l8L [l9|] . 



and the cosmic ray data may 
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the source because most of them are produced well outside our galaxy. With 
small but non-negligible probability the cosmic-ray protons can interact with 
EBL outside our galaxy, at distances (10-100) Mpc. The pions produced in 
the processes pj — > pre and pj — > nix + quickly decay and produce pho- 
tons and electrons with very high energies. These secondary photons can 
be observed by the Cherenkov telescopes, and they can create an image of 
the remote source even if all or most of the primary VHE photons produced 
by the AGN are filtered out by their interactions with EBL. Cosmic ray 



interactions with EBL should also produce neutrinos |21| . 

Let us consider the hypothesis that all the photons observed by the 
gamma-ray telescopes from distant AGN are, in fact, secondary photons 
from interactions of AGN produced cosmic rays with the background pho- 
tons. Production of cosmic rays is subject to model dependent assumptions, 
but it is generally believed that the output of a single AGN in cosmic rays 
with energies above 10 9 GeV is Fagn = (10 43 — lO 45 )^ 1 erg/s, where Eg 



is energy in units of 10 9 GeV 22]. The flux can be beamed with the beam- 
ing factor /beam ~ 10 — 10 3 , which increases the flux of protons from an 
AGN with a jet pointing in the direction of Earth (blazar). Some fraction 
/int ~ 10~ 3 of the protons can interact with EBL close to Earth and pro- 
duce pions, which generate photons and electrons as their decay products. 
These photons cascade down to energies below the pair production threshold 



and generate / mu i t ~ 10 3 photons below the threshold 23]. This series of 
interactions produces the flux 



F 1 ~ lO-^cm-V 1 



F AGN [E>10 7 GeV] \ 
10 46 erg/s j 



v (fb*m\ ( /int \ ( fmu\t \ ( \ Gpc \ 2 

which is can be large enough, and which is uncertain enough to be consis- 
tent with the fluxes detected by 7-ray telescopes. In fitting the 7-ray data 
we will use the spectral slope oc E~ 2 predicted by models of cosmic ray ac- 



celeration [22|, and we will choose the normalization based on the observed 
7-ray flux. We will not make assumptions regarding the luminosity of a given 
source in cosmic rays, and we will use the luminosity as a fitting; parameter, 



as long as it is consistent with the overall energetics of AGN |22J. We note, 
however, that harder intrinsic blazar spectra than those considered here can 
be produced by relativistic shock acceleration |24|. 
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Figure 1: The mean pion energy, and the corresponding mean photon energy, depend on 
the spectrum of EBL at a given red shift. Left panel shows n(e) for redshift z = (lower 
line) and z = 0.6 (upper line), according to Refs. 0, Q. The panel on the right shows the 
spectra of pions produced by interactions at z = 0.28 (the spectrum with a lower mean 
energy) and z = 0.014 (the spectrum with a higher mean) by protons emitted from the 
source at z = 0.444 (e.g., 3C66A) with a E~ 2 spectrum [22j. The units are arbitrary, and 
the upper curve is scaled down by a factor 10. 



We have calculated the spectrum of secondary 7-rays numerically using 
the EBL spectrum calculated in Refs. 0, Q (see Fig. 1) in the range z = 
— 0.6. The energy losses are due to production of pions in py interactions 
with the EBL photons. The interaction length A is given by [25] 



A- X =/ |^t/ a(s)(s-ml)dsde, (2) 

r>. wC 1—J J Sfn.i.'n. 
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where s min = ml + m 2 pl s max ~m 2 p + AeE p and e min ~ m 7r (m 7r + 2m p )/AE p . 
The photon density n(z, e) includes both the EBL and cosmic microwave 
background (CMB) photons, as shown in Fig. 1, and e max is the cutoff at 
about 160 eV. We calculated the optical depth, t(z), of the proton traveling 
from the source at z < 0.6 to z = in steps of ~1 Mpc. We assume that 
the protons emitted from AGN have spectrum oc E~ 2 at 10 6 — 10 9 GeV, 



consistent with some acceleration models [22|, and we use this spectrum 
as input for calculating numerically the rates of pion production at various 
values of the redshift. These pions, which carry about 12% of the proton 
energy, decay and initiate the electromagnetic cascades [23|. The showers 
produce both gamma rays and charged particles which should pair produce 
and undergo inverse Compton scattering (ICS) off the EBL and CMB until 
their products can propagate the distance to the earth. The photon flux is 
integrated numerically over the range of pion production redshifts. 
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Figure 2: Differential spectrum of secondary photons calculated numerically for 3C66A 
(z = 0.44), assuming the proton spectrum at the source oc E~ 2 , and normalized to fit 
MAGIC (blue squares) and VERITAS (red circles) data points. 



The predicted spectrum of 7-rays turns out to be similar for all the distant 
AGN. We have calculated the spectra for redshifts of 3C279, 1ES 1101-232, 
3C66A. The predicted spectra do not show significant differences because 
the two main factors that affect the spectrum of secondary photons, the 
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intrinsic spectrum of protons and the EBL, are similar for these sources. 
We show in Fig. 2 the spectrum of secondary VHE photons that should be 
observed from 3C66A (z = 0.44). The normalization is chosen to fit the 
data, and the spectral slope is in reasonable agreement with the observed 
spectrum of photons. Assuming the beaming factor /beam = 10 2 , the best- 
fit normalization of the flux implies the intrinsic luminosity in cosmic rays 
Fagn ~ 3 x 10 46 erg/s at energies above 10 7 GeV. Of course a higher fb ea m 
would imply a lower value of Fagn, and vice versa. The best fit to the 
spectral slope, dN/dE y = foE^, is T = —3.68 in agreement with MAGIC 
and VERITAS results. 

Our interpretation of the gamma-ray data is consistent with the obser- 
vations of point-like sources if the intergalactic magnetic fields (IGMF) are 
sufficiently small. Considerations based on Faraday rotation provide some 
upper limits on the magnetic fields in some directions, but there is no direct 
observational data for IGMF. Constrained simulations 171, ll9| predict a very 
inhomogeneous distribution of magnetic fields with voids, filaments, etc., but, 



as emphasized in Ref. [17J], the results of these models should be taken as 
upper limits on IGMF. The deflection angles depend on the distance to the 
source D, the correlation length l c , and the average value of the magnetic 
field B: 

It is possible that IGMF along the line of sight are as small as B ~ 10 -14 — 
10~ 13 G, in which case the image of a remote blazar in the secondary photons 
is point-like. Even if the line of sight crosses relatively narrow filaments or 
walls with a higher magnetic field, the angular displacement of the apparent 
image is not expected to be large, as long as such features are sufficiently 
distant. The magnetic fields in the voids can be well below B ~ 10 . 

The magnetic fields of the host galaxy do not affect the pointlike nature 
of the image. Even if the proton deflections in the host system are by large 
angles, A^host ~ 1, the enlargement of the image cannot exceed L^ ost /D <C 
0.1°, where Lh os t is the angular size of the host system. 

Of course, the cosmic ray detectors would not observe a point-like image 
at these energies because of the cosmic ray deflections caused by the galactic 
magnetic fields, which are much greater than IGMF. 

The photon-rich environments of M87 and 3C279 are not likely to allow 
the primary TeV photons to come out [3j. However, the cosmic rays are not 
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attenuated significantly by the same backgrounds. The observations of M87 
and 3C279 are consistent with the secondary photon interpretation. 



Multiwavelength observations of blazars [12| show the lack of correlation 
between the optical and VHE bands of variable sources. Furthermore, recent 
observations of M87 point to a significant mismatch of variability timescales 
of the likely emitting region observed in X-rays and in VHE photons I l3l| . If 
the VHE photons detected by HESS H, MAGIC 0, and VERITAS Q are 
secondary photons, then both the lack of time correlations and the differences 
in variability time scales can be explained by the delays in the arrival times 
of protons as compared to photons. These delays depend on the IGMF and 
on the proton energy, E p : 

,„4 I B \ 2 fW 7 GeV\ 2 ( D \ 2 I l c \ 

If IGMF are as low as 10~ 17 G, the energy-dependent time delays for the 
Markarian 501 flares would be of the order of a few minutes. Since Mrk 501 is 
relatively close [z = 0.031), we expect the low-energy signal to be dominated 
by primary gamma-rays, while the high-energy photons would be secondary 
photons, arriving with a delay. This is consistent with the MAGIC obser- 



vations of energy-dependent delays from a Mrk 501 flare [14(. If this is the 
case, no such delay should be seen for more distant sources, such as 3C66A 
and M87, where all the photons observed should be secondary photons. The 
assumption that IGMF are smaller than 10 -17 G is consistent with both 



the observational data and some of the models 281, while other models can 



accomodate larger IGMF [l8|, |29|, [l9|. 



The pair 3C66A and 3C66B provides an interesting case. 3C66A is be- 
lieved to be a blazar at distance 1.7 Gpc. The radio galaxy 3C66B appears 
to be only 6' away from the line of sight, at distance of 88 Mpc. VER- 
ITAS ji| has reported observations of VHE photons from 3C66A but not 
3C66B, while MAGIC has apparently observed photons from 3C66B 30). It 
is possible that the magnetic fields of the radio galaxy close to the line of 
sight could bend (and focus) the trajectories of the high-energy protons pro- 
duced by blazar 3C66A. Therefore, it is possible that MAGIC has observed a 
flare that originated in 3C66A as a burst of cosmic rays, which were focused 
by the radio galaxy, and which produced secondary photons pointing back 
to 3C66B, rather than 3C66A. 

Identification of detected photons with secondary photons from AGN al- 
lows one to use the gamma-ray data to study the EBL and the acceleration 
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mechanisms of cosmic rays. Although the observed spectra depend on the 
convolution of the two quantities, namely, the EBL and the cosmic ray spec- 
trum, the future data from similar sources at different red shifts may allow 
one to deconvolve the uncertainties in these two distributions. 

Future studies of gamma-ray data can test our interpretation of the pho- 
tons that appear to originate in remote AGN. On the one hand, it is possible 
that IGMF are high and EBL density is low [7[, in which case the primary 
gamma-rays from distant blazars are observed by gamma-ray telescopes. On 
the other hand, it is possible that IGMF are low, while EBL density is rel- 
atively high jl, 0, 0] , in which case the most distant AGN are seen almost 
entirely in secondary photons. If, indeed, the photons detected on Earth are 
secondary photons produced by the interactions of cosmic rays with EBL, 
this association will have a number of ramifications for understanding EBL, 
IGMF, and the mechanism of cosmic ray acceleration in AGN. 

The authors thank R. Ong, F. Stecker, and V. Vassiliev for very helpful 
comments, and M. Malkan for discussions of EBL and for providing some 
very useful data. This work was supported in part by DOE grant DE-FG03- 
91ER40662 and by the NASA ATFP grant NNX08AL48G. 
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